function [M_full,var_M_full,W_full,X_monthly,random_state] = data_moments(project_folder,which_moments)

Nu                                      = initialize_random_environment;
adoptionrates                           = importfile([project_folder '/input/' 'adoption_rates.csv'],2,7169);
X_monthly                               = adoptionrates.X_it;
X_monthly                               = reshape(X_monthly,[14 Nu.NumberOfDistricts]);
if nargin > 1 && strcmp(which_moments,'baseline')
    [M_full,var_M_full]                     = compute_data_moments(X_monthly,Nu);  
elseif nargin > 1 && strcmp(which_moments,'persist')
    [M_full,var_M_full]                     = compute_data_moments_IRF_cf_persist(X_monthly,Nu);  
else
    [M_full,var_M_full]                     = compute_data_moments_IRF(X_monthly,Nu);
end
W_full                                  = (( (numel(M_full)) )*( var_M_full ))^(-1); 
random_state                            = rng; 

end